Testing the QCD fragmentation mechanism on heavy quarkonium production at LHC 
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We calculate the fragmentation function for charm quark into J /if) at the QCD next-to-leading- 
order (NLO) and find that the produced J/ip is of larger momentum fraction than it is at the 
leading-order. Based on the fragmentation function and partonic processes calculated at the NLO, 
the transverse momentum distribution on J/ip hadroproduction associated with a charm c (or c) 
jet are predicted. We find that the distribution is enhanced by a factor of 2.0~3.3 at the NLO 
as pt increased from 10 GeV to 100 GeV and it is measurable at the LHC with charm tagger. 
The measurement at the LHC will supply a first chance to directly test the QCD fragmentation 
mechanism on heavy quarkonium production where the fragmentation function is calculable in 
perturbative QCD. It is also applied to J/ip (T) production in the decay of Z° (top quark). 
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Quantum Chromodynamics (QCD) is a successful the- 
ory to describe strong interaction, but its fundamen- 
tal ingredients, the quarks and gluons, are not observed 
freely and must hadronize eventually. The fact makes it 
impossible to calculate any processes involving detected 
hadrons in the final or initial states directly. According 
to the QCD factorization theorem (see Ref. [l[ and ref- 
erences therein), in some kinematical regions, the domi- 
nant contribution to the cross section can be decomposed 
into three parts: the partonic part, the part of parton 
fragmentation into the produced hadron and the part of 
the parton distributions in the initial hadrons. The par- 
tonic part can be calculated perturbatively because of the 
asymptotic freedom of QCD, while all the long distance 
physics of the hadrons is put into the parton fragmenta- 
tion functions and the parton distributions. Therefore, 
fragmentation functions is one of the most important in- 
gredients to understand QCD or to make predictions for 
experimental measurements. It is hard to study frag- 
mentation functions directly from QCD because of their 
non-perturbative nature. 

For light hadrons, the fragmentation functions are ex- 
tracted from global data fits. Recently the transverse- 
momentum pt distribution of inclusive light-charged- 
particle production measured by the CDF shows sig- 
nificantly exceed on the theoretical prediction based on 
these fragmentation functions when p t > 80 GeV @ • It 
potentially challenges our understanding of QCD fac- 
torization theorem. However, for heavy quarkonium, 
the non-relativistic QCD (NRQCD) factorization formal- 
ism 0] can be used to factorize the fragmentation func- 
tions for quarkonium into NRQCD matrix elements and 
short-distance factors, which are calculable in perturba- 
tion theory and have been studied in many works [H-[6[ 
at QCD leading-order (LO). There are also QCD next- 
to-leading order (NLO) studies on the color-octet 3 Si 
gluon fragmentation function for heavy quarkonium Q, 
and the study on relativistic corrections for the frag- 
mentation functions @. The QCD factorization on 



heavy quarkonium production is investigated by many 
authors (^.Therefore, we have more prediction power to 
test QCD fragmentation mechanism on heavy quarko- 
nium production. But until now, there is no available 
experimental measurement to test it. 

In recent years, there is a huge data collection at col- 
liders. Based on that, many J/ip production processes 
were observed [To[ in the past. Therefore it supplies a 
very important chance to perform systematical study on 
J/ip production both theoretically and experimentally. 
The large discrepancies for exclusive J/ip productions at 
the B factories have been studied and resolved by in- 
troducing higher order corrections [ll[. The discrep- 
ancies for inclusive J/ip production at the B factories 
have also been studied [lj, EH and the results including 
NLO QCD corrections can nearly explain the experimen- 
tal data. Higher order corrections for J/ip production at 
hadron colliders were also investigated [Tj]. Although 
large improvement have been achieved in theoretical pre- 
dictions, the experimental data are still unable to fully 
understand, especially for the polarization of J/ip. 

Higher order contributions have shown their impor- 
tance and the dominant production mechanism for heavy 
quarkonium is fragmentation at large transverse momen- 
tum region. To achieve reasonable theoretical predic- 
tions for J/ip production at large pt, it is important to 
study higher order contribution to fragmentation func- 
tions for heavy quarkonium states. Moreover, the most 
important question here is: is there any chance to di- 
rectly test the QCD fragmentation mechanism on heavy 
quarkonium production at the LHC? In this letter, we 
calculate the fragmentation functions for charm quark 
into J/ip at QCD NLO. Based on the fragmentation func- 
tion and partonic processes calculated at the NLO, the 
transverse momentum distribution on J/ip production as- 
sociated with a charm c (or c) jet are predicted. We find 
that the distribution is enhanced by a factor of 2.0^3.3 
at the NLO as p t increased from 10 GeV to 100 GeV and 
it is measurable at the LHC with charm tagger [l5[ . 
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Here we study J/ip production in e + e~ — > J/ijj + X 
with a very high cm. energy y/s, which should be dom- 
inated by fragmentation mechanism. The differential 
cross section for J J ip with momentum p is 

da[e + e--> J/^(p) + X) (1) 
= Jdzda{e + e~^i(p/z)+X,Li F }D z ^j/ 4 ,(z,n F ). 




7 r" X 






;qrt = 




sqrt{s)--0OGeW 




5qrt{s)=3G0GcV 









Where the factorization scale [ip is introduced to main- 
tain this factorizcd form and the dependence on the ar- 
bitrary scale fiF cancels between the two factors. And i 
could represent all quarks and gluon, but we only con- 
sider i = c (i.e. e + e" — > J/'tp + cc channel) in the theoret- 
ical calculation for both side of Eq.(l) to extract charm 
quark c to J/ip fragmentation function. Here we choose 
1-iF = 3to c to avoid large logarithms of \xf jm c in the frag- 
mentation function D c ^_,.j^(z, Hereafter D(z) is 
used to represent D(z,3m c ). Thus we have 
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up to NLO in a s with z = Ej/^,/E c . As mentioned 
in Refs. [3, [H, the initial fragmentation function can be 
calculated perturbatively as a series in a s (2m c ) and ex- 
tracted from Eq. ((2|) order by order. 

At leading-order (LO) in a s , there are E c = E s = y/s/2 
and z = 2Ej/^,/y/s, and Eq. ([2]) is simplified into 



e+e-^J/tjj+X 



LO 



dE 



j/i> 
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Thus the fragmentation function is extracted as: 
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Under the limitation m c /y/s -> 0, a* = (647ra 2 )/(9s 3 / 2 ) 
is obtained, and by using Eq.(4) in Rcf. [13| we obtain 
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which is exactly the same as the one obtain from Z decay 
in Rcf. |[ . 

In the fix order calculation at NLO in a Si Eq. ([2]) be- 
comes 
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FIG. 1: Behavior of fi{z) and fc{z) with fiF = 3m c . 

Then the NLO fragmentation function is expressed as 

D^?/4z) = h{z)-h{z) (7) 
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fi is the energy distribution of J/ ipec production at NLO, 
which have been achieved in our previous work. And /a 
needs the LO fragmentation function in Eq. (JSJ) - It is hard 
to obtain an analytic result here and we have to do it nu- 
merically. We choose m c = 1.5 GeV, a s {2m c ) = 0.26 
and |i? s (0)| 2 = 0.944 GeV, and the renormalization scale 
/jlr = 2m c . The behaviors of fi(z) and f2(z) shown in 
Fig. Q] strongly depend on y/s. In Fig. [2^, the fragmen- 
tation functions extracted numerically at the LO and 
NLO show very good limitation as y/s increasing from 
30 GeV to 1000 GeV and even the difference between 
the y/s~ = 30 GeV and y/s~ = 1000 GeV is quite small. 
It means that the fragmentation mechanism can describe 
the theoretical result even when the cm. energy is as low 
as 30 GeV. For the NLO results, there exists unphysical 
range of negative possibility arising from fix-order per- 
turbative calculation. The result with y/s = 1000 GeV 
is a bit unstable even in our quadruple precision FOR- 
TRAN calculation because of large numerical cancella- 
tion. Therefore we choose the result at y/s = 300 GeV 
to approximate the final D^^^(z). It is clearly shown 
in Fig. [3] that J/tp from the fragmentation is of larger 
momentum fraction at the NLO than that at the LO. 



On the other hand, D^^ Uh {z) can be expressed as 



n NLO ( 
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a(z) + /3 ln — — 
2m c 



where (3q is the one-loop coefficient of QCD beta func- 
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tion. A 9th-polynomial fitting gives a(z) = qz' with 

Co = -4.14253418603 x 10\ci = +1.08551074161 x 10 2 , 
c 2 = +1.62953354180 x 10 3 ,c 3 = -1.77039359042 x 10 4 , 
c 4 = +7.69974017791 x 10 4 ,c 5 = -1.82373533780 x 10 5 , 
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And the two corresponding NLO results are 



FIG. 2: The fragmentation functions at LO and NLO (In left 
figure, jj,F—3m c is chosen and lower curves are LO ones). 



c 6 = +2.52868747876 x 10 5 ,c 7 = -2.03305581610 x 10 5 , 
c 8 = +8.66677519749 x 10 4 ,c 9 = -1.47927072151 x 10 4 . 
This is the NLO fragmentation function of charm into 
J /if) at the initial factorization scale [if = 3m c . The 
fragmentation function at other factorization scale can 
be obtained with this initial one by solving the Altarelli- 
Parisi evolution equation [l6| . In Fig.[2p the evolutions of 
the LO and NLO fragmentation functions are presented. 

It is easy to do applications with this new fragmen- 
tation function. Firstly we study the decay of Z° into 
J ftp in Z° — > J/ip + cc. With the additional parameter 
a = 1/128 , the decay width at LO is 



^LO 

J/4>+x 



2T 



LO 
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dzD 



LO 



(z) = 129 KeV. 



At NLO, there are two ways to calculate the decay width. 
One is described by Eq. where the higher order term 
is neglected 



yNLO _ 9J ,LO I f l 7 r,NLO / \ 
1 J/ip+X — ZL c+X I azU c^J/il>\ z ) 



c+X 
dT NLO- 



dE c dz J+ x (z) = 136 KeV. 

The other one is to include the higher order term as 



dT 



NLO 

c+X T-.NLO 



D^Y/^z) = 141 KeV. 



Both the LO and NLO results are consistent with the 
fully NLO QCD calculation 0, which gives 120 KeV at 
LO and 136 KeV at NLO with same parameters. The 
differences come from the fact that the limitation is not 
so well as the mass of Z° is not large enough to be treated 
as infinity. Secondly we apply the fragmentation function 
to b quark case by substituting 

rn c ++m b , n f = 4 -H- n f = 5, i? s 7/!/> (0) ++ Rj(0). 

For the top quark decay, t — > T + W + + b, we have 

T^t+x = T t°b+x J dzD£° r (z) = 30.9 KeV. (9) 



r NLO 

1 t^r+x 



T^+x I dzD^°(z) 



'b->r 



dT?A°* 



dE b dz — l /^ +x D%° r (z) = 40.0 KeV, 



dE b 



^NLO+ 
t^T+X 



= / dEhdz- 



dY 



NLO* 

t^b+X r^NLO 



dE h 



Db-Xr 0) = 39.7 KeV. 



Here we choose the same parameters as those used in 
Ref. [l8|]. One should notice that the LO wave function 
at the origin is used in Eq. © , while in other two results 
the NLO one is used. The corresponding LO and NLO 
results given by Ref. [HI are 26.8 and 52.3 KeV. 

The most important application is to study the produc- 
tion of J/tp + cc + X at the LHC. Because the fragmen- 
tation function Dg^j^ is suppressed comparing with 
D c (c)^-j/if>i we only consider the contribution from the 
charm quark fragmentation. The cross section is 

<j[pp — > J/ipcc + X] 
= ^2 J dxidx 2 dzf i / p (xi^f)f j / p (x2,Hf) (10) 

x da[ij -+ cc + X,n f ,n r , £iF]-D c (e)->j/i/> (z, fi F ) , 

where f(x, [if) is the parton distribution function, [if, [i r 
and [if are the factorization, rcnormalization and frag- 
mentation scales respectively. And a represents the cross 
sections of partonic process. The LO and NLO fragmen- 
tation functions are used to calculate the final LO and 
NLO results respectively. In the calculation, there are 
two subprocesses gg{qq) — > cc at the LO, and three real 
corrections gg — > cc + g, gq{q) — >■ cc + q(q), qq — > cc + g 
and the virtual corrections to gg(qq) — > cc at the NLO. 
The default choice of charm quark mass is m c = 1.5GeV 
and the three scales are set as [if = [i r = [ip = [i 

with the default choice [i = [iq = J + m*. We 
choose m c = 1.4 GeV, 1.6 GeV and [i = [1q/2, 2[i for 
uncertainty estimation. The Cteq6Ll and Cteq6M [3] 
are used in the LO and NLO calculations respectively, 
with the corresponding a s running formula being used. 
The fragmentation function are evolved from 3m c to [If 
by solving the Altarclli-Parisi equation numerically. The 
wave function at the origin of Jhp is extracted from its 



l_7f at the NLO level and 
< 2.4. 



lcptonic decay as in reference 
the rapidity cut for J/tjj is \yj/$ 

The predictions for transverse momentum distribu- 
tions of J/ip are shown in Fig. [3] with t/s = 7 TeV and 
^/s = 14 TeV at the LHC. The complete calculation of 
process pp — > J/ip + cc + X at the LO are presented 
for comparison, which was studied in reference [2(|. It 
will be dominated over as pt > 40 GeV by the LO frag- 
mentation result and as pt > 16 GeV by the NLO frag- 
mentation one, and is about 5 times smaller than the 
NLO fragmentation one in large pt region. It is believed 
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(a) P, (GeV) 10 ™ °° ° (b) P, (GeV) 



FIG. 3: The theoretical prediction on pt distribution of J/tp 
production associated with a charm c (c) jet at the LHC. 



that the fragmentation mechanism should give a better 
description on the J/tp hadroproduction associated with 
a charm c (c) jet at large pt region than the complete 
calculation result. 

We extracted the fragmentation function for charm 
into J/tp at the NLO in a s , D^°( z , 3m c ), from J/ip 
production at e + e~ annihilation. The fragmentation 
function for bottom into T, D^y(z,3mb), is also ob- 
tained. The Altarelli-Parisi evolution of the fragmenta- 
tion function is performed to obtain its value at other 
fragmentation scale. We have applied them to J/tp pro- 



duction in Z° decay, and T production in top quark de- 
cay. The most important application is to predict the 
transverse momentum distribution of J / tp production as- 
sociated with a charm c (or c) jet at the LHC. We find 
that the distribution is enhanced by a factor of 2.0~3.3 
at the NLO as pt increased from 10 GeV to 100 GeV and 
it is measurable at the LHC with charm tagger. In the 
measurement, it need to identify a J/tp and a charm jet 
without the J/tp in the jet and there will be 20 events 
to be found if we optimistically assume a 60% charm 
tagging efficiency, both e and \i decay channel of J/ip 
being detected and 10/& -1 of the integrated luminosity 
at a 14 TeV machine. Of course, this measurement is 
a very big challenge to the experimental technique on 
charm tagger. However, the measurement at the LHC 
will supply a first chance to directly test the QCD frag- 
mentation mechanism on heavy quarkonium production 
where the fragmentation function is calculable in pertur- 
bativc QCD. 
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